Phonon spectrum and soft-mode behavior of MgCNi 3 
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Temperature dependent inelastic neutron-scattering measurements of the generalized phonon 
density-of-states for superconducting MgCNi3, T c = 8K, give evidence for a soft-mode behavior 
of low-frequency Ni phonon modes. Results are compared with ab initio density functional calcula- 
tions which suggest an incipient lattice instability of the stoichiometric compound with respect to 
Ni vibrations orthogonal to the Ni-C bond direction. 

PACS numbers: 74.25.Kc,63.20.Kr,78.70.Nx,74.70.Dd 
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The recent discovery of superconductivity near 8 K 
in MgCNi^i was a surprise since due to the high Ni 
content a ferromagnetic ground state would have been 
expected rather than a superconducting one. Further- 
more it is unusual that for this perovskite compound 
the heavy Ni atoms occupy the corners in Ni 6 C octa- 
hedra. There are speculations that magnetic interac- 
tions might promote superconductivity^, however, mea- 
surements of the specific heat come up with the conclu- 
sion that MgCNi3 is a medium to strong coupling BCS 
superconductor^^. This compound bears similarity to 
superconducting nickel borocarbides where for the Y and 
Lu compounds an acoustic phonon branch exhibits a pro- 
nounced softening at lower temperature s - ' a feature 
which is supposed to be connected to superconductivity 
in these compounds. Hints for unusual lattice dynamical 
properties of MgCNi3 have been inferred from a recent 
x-ray absorption study, which indicated deviations of the 
local atomic structure from the ideal perovskite lattice at 
temperatures below 70 and from a very recent linear 
muffin-tin orbitals (LMTO) calculation of the harmonic 
phonons of MgCNigS, which found a dynamical instabil- 
ity of the stoichiometric compound. In this paper we will 
present inelastic neutron scattering (INS) measurements 
of the generalized phonon density-of-states (GDOS) as a 
function of temperature which show a softening of low- 
frequency Ni modes. We will compare our experimental 
results to ab initio density functional calculations. 

For the synthesis of MgCNi3 samples mixed powders of 
the components were pressed into pellets, wrapped into 
a tantal foil, sealed in a quartz tube and heated for 2 h at 
900° C. A difficulty arises from losses of volatile Mg and 
C at the synthesis temperature. With respect to subse- 
quent x-ray analysis we obtained the best samples for a 
starting stoichiomctry of Mg1.14C1.4Ni3 where an impu- 
rity peak due to unreacted graphite was smallest. The T c 
of our samples was 7.5 K. It could be slightly increased for 
the starting concentration C1.5 however at the expense 
of a somewhat larger contribution of unreacted graphite. 
Since the exact stoichiometry might be important for a 
soft-mode behavior which is found in the present inves- 
tigation we show in Fig. ^ a x-ray diffraction pattern of 
our sample. 
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FIG. 1: Observed (crosses) and calculated (solid line) x-ray 
diffraction pattern for the investigated sample. A small im- 
purity peak at 26.5° is due to unreacted graphite. The re- 
finement (R = 5%) yields a composition of MgC.9eNi3 for the 
cubic perovskite structure, space group Pm3m, with lattice 
constant a=3. 8085(1) A and with positions of the atoms at: 
Mg la (0,0,0), C lb (0.5,0.5,0.5), and Ni 3c (0,0.5,0.5). 



Our INS measurements were performed at the HFR in 
Grenoble, France, on the IN6 and IN4 time-of-flight spec- 
trometers. On IN6 a high chopper speed of 12060 rpm 
and focusing in the inelastic region were used to im- 
prove the resolution. Due to the low incident energy of 
4.75 meV at this spectrometer which is connected to the 
cold source we profit from a very good resolution for low- 
energy excitations. The necessity to work in energy gain 
of the scattered neutrons sets limits on the possibility 
to work at lower temperatures. For example, for a sam- 
ple temperature of 50 K we could follow excitations up 
to ~ 30 meV. The spectrometer is dedicated to inelastic 
measurements and a standard result is the GDOS which 
is calculated from the recorded intensities over a scatter- 
ing region from 14° to 114°. For the data evaluation we 
have applied multi-phonon corrections in a self-consistent 
procedure^. The GDOS implies a weighting of vibra- 
tional modes by a/m (scattering cross-section over the 
mass; we used values of 0.46, 0.179, and 0.31 barn/amu 
for C, Mg, and Ni, respectively). From a plot of the elas- 
tic intensities over momentum transfer we do not obtain 
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FIG. 2: The experimental generalized phonon density-of- 
states for MgCNi3. A significant soft-mode behavior is found 
for low- frequency Ni modes (insert). 



any evidence for a structural phase transition within the 
investigated temperature region in agreement with previ- 
ous results from structural invest igation o 1 ' 1 1 . However, 
due to the low incident energy this analysis is limited 
to Q <2.2A _1 . To study changes in the dynamics in 
a larger temperature region, we have performed supple- 
mentary measurements with an incident neutron energy 
of 35 meV in the down-scattering mode over a wide range 
of temperatures on the instrument IN4. These results 
will be presented in the form of the dynamical suscepti- 
bility x"( w ) which is much closer related to the scattering 
function S(Q, w) than the GDOSiS. 

Due to the large differences in mass the phonon spec- 
trum of MgCNi3 (Fig. |2J) decomposes into two well sep- 
arated parts, a low-frequency region with predominantly 
Ni and Mg modes and vibrations of the light C atom 
around 80 meV. From the known a/m values and an anal- 
ysis of the areas in the GDOS shown in Fig. [5] we find for 
the latter peak an almost perfect agreement between the 
theoretical and the experimental spectral weight imply- 
ing contributions from three closely neighbored optical 
phonon branches. The weight of the sharp and well sep- 
arated peak at 35meV corresponds well to the weight 
expected for three optical Mg modes. It is increased by 
some 20% which suggests a small hybridization with the 
low-frequency Ni modes. The region below 33 meV con- 
tains exclusively Ni contributions. It is well structured 
and exhibits a very low frequency maximum at 12meV 
and a strong main peak around 16 meV. 

The 12 zone center optic phonons decompose according 
to 3F lu + F2 U with no Raman active modes. Guided by 
investigations of the related system of quaternary boro- 
carbides where a soft-mode behavior for superconducting 
RENi2B2C (RE = Y, Lu) has been found we have mea- 
sured the GDOS of MgCNi3 for a series of temperatures 
between 400 K and 50 K. No changes are found for the 
peaks at 80 meV and 35 meV which were attributed to 
C and Mg vibrations (only a small sharpening of the 
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FIG. 3: The experimental dynamical susceptibility for 
MgCNi 3 (see text). 



35 meV peak is registered) . A pronounced softening, 
however, occurs for the very low-frequency part of the 
Ni vibrations in a region which is dominated by acoustic 
phonon branches. The insert in Fig. [5] shows a magni- 
fication. The strongest effect is observed between 200 K 
and 100 K. Only less pronounced changes occur between 
100 K and 50 K which is the lowest temperature in our 
energy-gain measurements. 

Changes in the dynamics show up clearer in the tem- 
perature evolution of x"(w) in Fig. [21 where a sum of 
intensities recorded on IN4 for different scattering angles 
is shown. The incident neutron energy was chosen such 
that the low-energy region below 20meV could still be 
studied with good resolution. The strong increase in in- 
tensity below 3 meV is due to elastic scattering contami- 
nations. Whereas the frequency of the prominent Ni peak 
remains almost fixed at 16meV we can observe the evo- 
lution of the soft-mode intensity. In agreement with the 
GDOS the strongest changes are observed around 100 K 
and only minor changes are recorded between 300 K and 
200 K and between 30 K and 2K. Additionally to the 
GDOS in Fig. [21 it can be seen that some intensity is 
shifted into the low-frequency region below 5 meV. While 
the increase in intensity below 10 meV is in general agree- 
ment with the mode softening a more detailed analysis of 
its location within the Brillouin zone (BZ) is not possible 
for this polycrystalline sample. 

For an analysis of our data we have performed first 
principals density functional calculations using a mixed 
basis pseudopotential method. For Mg a well tested pseu- 
dopotential of Troullier-Martins type has been usediSiH, 
whereas for C and Ni Bachelet-Hamann-Schliiter type 
pseudopotentials were constructed^. The fairly deep d 
potential of Ni and p potential of C can be efficiently 
dealt with by the mixed-basis formalism, where the va- 
lence states are constructed from a combination of local- 
ized s and p functions at C sites and localized d functions 
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FIG. 4: Calculated GDOS spectra for stoichiometric MgCNi 3 . 
The proper a/m values have been applied for a comparison to 
the experimental data (300 K spectrum). Contributions from 
an unstable phonon branch (see text) have been ignored. 



at Ni sites, supplemented by plane waves up to a kinetic 
energy of 24 Ry. The local density approximation with 
the Hedin-Lundqvist parametrization of the exchange- 
correlation functional^ has been used. Our calculated 
electronic band structure is in good agreement with pre- 
vious density functional studiea 2 i 16 i 17 i 18 i 19 i 20 i 21 i 22 . The 
most prominent feature of the electronic density-of-states 
is the existence of a strong and sharp van Hove peak 
less than 100 meV below Ef ■ It arises from antibonding 
bands related to the Ni d orbitals, which are very narrow 
due the linear Ni coordination by G™. 

Phonon dispersion curves for stoichiometric MgCNi3 
were calculated within the mixed-basis perturbation 
approach^. Complete spectra are obtained from a 
Fourier interpolation of dynamical matrices calculated on 
a cubic (4x4x4) q-point mesh. Results for the optimized 
lattice constant a=3.76A are depicted in Figs. 01 and [5] 
The theoretical GDOS (Fig.^J reproduces many features 
seen in the experimental spectra, including the small 
maximum near 30meV and the shoulder of the acoustic 
spectrum at 20meV. The acoustic modes are predomi- 
nantly of Ni character with only a small hybridization of 
C modes. The sharp peak near 35 meV arises from three 
almost dispersionless Mg branches with a small admix- 
ture of Ni vibrations. 

The most prominent feature of the calculated disper- 
sion curves (Fig.[3J) is, however, the occurrence of unsta- 
ble modes suggesting that the stoichiometric compound 
is dynamically unstable in the harmonic approximation. 
This instability appears over a large area in reciprocal 
space and is most pronounced near the zone boundary 
points X and M (with coordinates (100) and (110) in 
units of 7r/a, respectively). Inspection of the displace- 
ment vectors reveals that all unstable modes involve Ni 
vibrations orthogonal to the Ni-C bond directions. This 
behavior is most evident at the M point where the unsta- 
ble mode represents a pure vibration of the Ni atoms with 




FIG. 5: Calculated phonon dispersion curves for stoichiomet- 
ric MgCNi3 for selected symmetry directions. Black circles 
denote the results from linear-response calculations. Lines are 
obtained by Fourier interpolation. Negative frequencies indi- 
cate unstable modes. At the M point, the unstable phonon 
corresponds to shear elongations of Ni atoms in the Ni-C 
planes as indicated in the insert. 



all other atoms at rest. Its elongation pattern, shown in 
the insert of Fig.O induces shear distortions in the planar 
Ni-C sublattice. This extended phonon instability may 
be the origin of the apparent downshift of the acoustic 
part below 20 meV in the calculated spectrum by 3 meV 
with respect to the experimental one (see Fig.0J. 

The presence of the van Hove peak close to Ep re- 
quires care in performing the BZ summations in the 
calculations. The presented results were obtained with 
cubic (16x16x16) Monkhorst Pack k-point meshes and 
a Gaussian broadening of 0.1 eV. Convergence studies 
showed that this is sufficient for most of the phonon 
modes to be converged within 0.5 meV, but the low- 
frequency Ni modes are more sensitive to the BZ 
sampling. Frozen-phonon studies of the unstable M- 
point mode using a variety of k meshes and Gaus- 
sian smearings have confirmed that this instability is 
not the result of an insufficient BZ sampling. Fur- 
thermore, a similar frozen-phonon calculation performed 
with the projector-augmented plane wave approach as 
implemented in the Vienna ab-initio simulation pack- 
age (VASP-PAW) 24 i 25 i 26 i 27 reproduced the instability, 
excluding the possibility that it is an artefact of the pseu- 
dopotential approximation. 

An instability of Ni modes has been also found in 
a very recent ab initio investigation of the lattice dy- 
namics and electron-phonon coupling (EPC) of MgCNi3 
using the LMTO method 9 . In contrast to our calcula- 
tions, however, the instability did not appear at high- 
symmetry points at the BZ boundary. While both meth- 
ods agree with respect to the frequency range for the 
Ni and C modes, the LMTO calculations predict signifi- 
cantly higher frequencies for the Mg modes (at «45 meV) 
at variance with our experimental and theoretical results. 



4 



The instability of the harmonic lattice dynamics of sto- 
ichiometric MgCNi 3 raises the question of the true low- 
temperature structure of this superconductor. Some in- 
sight can be gained from frozen-phonon studies for the 
unstable mode. For the shear mode at M we found a 
very shallow double-well potential, with a minimum of 
less than 1 meV at a displacement of only 0.025 A. In the 
previous LMTO study, double-well potentials of similar 
magnitude have been extracted from frozen-phonon cal- 
culations corresponding to lower-symmetry points 9 ". The 
shallowness of the double well suggests that the instabil- 
ity is not strong enough to induce a long-ranged struc- 
tural distortion, but that the high-symmetry cubic struc- 
ture is stabilized dynamically. The presence of strong an- 
harmonic motions of atoms in the vicinity of a structural 
instability could be of relevance for superconductivity, 
as it may significantly enhance the EPC strength. For 
MgCNi3, a contribution of 37% coming from the anhar- 
monic modes to the total EPC constant has been esti- 
mated on the basis of the LMTO results 9 . A complete 
calculation of the EPC strength, however, is a formidable 
task as it requires to take into account the full anhar- 
monic lattice dynamics as well as anharmonic corrections 
to the EPC. This is further complicated by the fact that 
the local distortions are expected to strongly modify the 
van Hove singularity and thus the electronic structure in 
the vicinity of the Fermi energy^. 

The picture of a dynamically stabilized cubic struc- 
ture is consistent with published structural investigations 



where no indications of any phase transition have been 
found, and with the temperature dependence of the elas- 
tic scattering intensity observed in our IN6 spectra. In a 
recent x-ray absorption fine structure measurement, an 
unusually broad Ni-Ni pair-correlation function has been 
observed below 70 KS. It was interpreted as the result of 
local distortions of the Ni octahedra predominantly per- 
pendicular to the Ni-C bond with displacements of less 
than 0.05 A, in agreement with the frozen-phonon results. 

In conclusion, studies of the GDOS have provided a 
clear insight in the lattice vibrations of superconducting 
MgCNi3 which can be looked at as a 3D analogue of the 
nickel borocarbides. For both kinds of compounds a re- 
markable soft-mode behavior in the low-frequency region 
(~ 8 meV) is observed. For the superconducting borocar- 
bides where single crystals could be grown it was shown 
that the mode softening is restricted to a narrow region in 
reciprocal space and is sensitive to the formation of the 
superconducting gapi In contrast, the present results 
suggest that the softening phenomena in MgCNi 3 is not 
directly linked to the superconducting state, but has its 
origin in an incipient dynamical instability of specific Ni 
vibrations which occurs widespread in reciprocal space. 
It cannot be excluded, however, that lattice defects play 
an important role in stabilizing the real structure. Fur- 
ther studies are required to clarify the question whether 
the apparently large stability range of the MgCi-^Nia 
phase is connected to the observed dynamical peculiari- 
ties. 
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